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Abstract 

This  work  presents  measurements  of  neutral  xenon  velocity  in  the  plume  of  the  Busek  Company  BHT-200- 
X3  200  W  Hall  thruster  at  a  single  250  V  anode  potential  operating  condition.  The  xenon  velocities  were  measured 
using  laser  induced  fluorescence  of  the  6s’[  l/2]j  -6p’[3/2]2  excited  state  transition  at  834.68  nm.  Velocities  were 
interrogated  near  the  exit  plane  of  the  anode  discharge  and  in  the  vicinity  of  the  cathode.  Both  axial  and  radial  veloci¬ 
ties  were  measured  simultaneously  using  phase  sensitive  detection.  The  exit  plane  axial  velocity  varies  from  350-50 
m/s  near  the  outer  acceleration  channel  radius  to  approximately  150  m/s  at  the  inner  radius.  The  radial  velocity  at  the 
outer  insulator  is  approximately  zero,  but  increases  to  approximately  100  m/s  directed  away  from  the  protruding  cen¬ 
tral  magnetic  pole  at  the  inner  insulator.  This  is  consistent  with  the  creation  of  neutrals  at  the  inner  radius  by  ions  neu¬ 
tralized  following  a  wall  collision  with  the  protruding  central  magnetic  pole.  Neutral  velocities  downstream  of  the  exit 
plane  also  show  a  radial  velocity  away  from  the  central  magnetic  pole.  However,  this  behavior  may  be  due  to  the 
influence  of  the  cathode  plume.  Velocity  measurements  in  the  cathode  plume  show  that  xenon  velocities  reach  in 
excess  of  610  m/s.  Since  the  cathode  flow  represents  a  substantial  portion  of  the  thruster  flow  and  is  only  weakly  ion¬ 
ized,  the  cathode  plume  may  influence  the  neutral  velocities  nearer  the  anode.  Evidence  of  interaction  between  the 
two  flows  takes  the  form  of  widened  fluorescence  line  profiles  which  may  indicate  the  collisional  mixing  disparate 
neutral  streams  or  possibly  several  gas  dynamic  shocks. 

Introduction 

The  goal  of  this  study  is  to  measure  the  neutral 
xenon  (Xe  I)  velocity  field  of  a  200  W  Hall  thruster 
using  excited  state  atomic  xenon  laser  induced  fluores¬ 
cence  (LIF).  Unlike  ions,  atoms  can  not  be  characterized 
by  electrostatic  probes.  Often,  LIF  is  the  preferred  diag¬ 
nostic  for  the  investigation  of  various  plasma  species 
since  it  is  a  non-intrusive  measurement.  Unlike  electro¬ 
static  probes,  LIF  will  not  perturb  the  local  plasma,  or 
thruster  operation.  In  this  particular  case,  LIF  appears  to 
be  the  only  practical  measurement  method. 

Hall  thrusters  produce  thrust  by  the  accelera¬ 
tion  of  ions  through  an  electric  field.  Xenon  is  used  as  a 
propellant  due  to  its  high  atomic  mass,  inert  chemistry, 
and  low  ionization  energy.  Xenon  ions  will  exit  Hall 
thrusters  at  ~15  km/s.  A  fraction  of  the  xenon  propellant 
is  not  ionized,  and  is  not  accelerated.  A  well  known  loss 
mechanism,  these  neutrals  probably  represent  less  than 
10%  of  the  propellant  flow.  However,  their  density  is 
comparable  to  ionic  flow  due  to  the  wide  difference  in 
velocities. 
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Understanding  the  behavior  of  the  neutral  flow 
is  important  for  the  complete  characterization  of  Hall 
thrusters.  For  example,  high  neutral  densities  influence 
electron  transport  by  providing  collision  partners  for 
low  energy  electrons  seeking  to  cross  magnetic  field 
lines.  Neutrals  are  also  a  source  of  charge  exchange  col¬ 
lision  partners  with  the  ion  stream.  Charge  exchange 
may  be  responsible  for  the  creation  of  highly  divergent 
low  energy  ions  which  influence  Hall  thruster  spacecraft 
integration.  In  addition,  how  the  neutrals  flowing 
though  the  thruster  interact  with  the  background  neu¬ 
trals  in  ground  test  facilities  is  a  growing  concern. 
Understanding  the  differences  between  ground  test  and 
on-orbit  thruster  operation  is  an  emerging  area  of  study. 

Complicating  neutral  flow  measurements  from 
the  anode  discharge  plume  is  a  high  density  jet  of  neu¬ 
trals  from  the  hollow  cathode  neutralizer  which  repre¬ 
sents  approximately  10%  of  the  total  thruster  propellant 
flow.  Unlike  the  anode  flow,  the  hollow  cathode  plume 
is  only  weakly  ionized.  The  cathode  produces  a  high 
density  jet  of  neutrals  which  may  also  influence  the  flow 
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field  in  the  vicinity  of  the  anode.  This  flow  of  neutrals 
may  affect  local  electron  conductivity  and  hence  plume 
neutralization.  As  a  result,  the  cathode  plume  may  also 
influence  insulator  erosion  patterns,  particularly  in  the 
case  of  very  high  power  Hall  thrusters. 

Xenon  Spectroscopy 

The  nine  stable  isotopes  of  xenon  each  have  a 
slight  difference  in  their  electron  transition  energies  due 
to  their  differences  in  mass.  The  odd  mass  isotopes  are 
further  spin  split  due  to  nuclear  magnetic  dipole  and 
electric  quadrapole  moments.  Nuclei  which  have  an  odd 
number  of  protons  and/or  an  odd  number  of  neutrons 
possess  an  intrinsic  nuclear  spin  Ih/2n  ,  where  I  is  inte¬ 
gral  or  half-integral  depending  on  whether  the  atomic 
mass  is  even  or  odd,  respectively.  For  nuclei  with  non¬ 
zero  nuclear  spin  (angular  momentum),  there  exists  an 
interaction  of  the  nucleus  with  the  electron  shell.  This 
interaction  leads  to  the  splitting  of  levels  with  angular 
momentum  J  into  a  number  of  components,  each  corre¬ 
sponding  to  a  specific  value  of  the  total  angular  momen¬ 
tum  F .  The  details  of  how  hyperfine  spin  splitting 
affects  the  spectra  of  xenon  are  described  in  more  detail 
elsewhere  [1]. 

For  the  results  reported  here,  the 
6s’[l/2\®  -6p’[3/2]2  electronic  transition  of  Xe  I  at 
834.68  nm  is  probed.  The  isotopic  and  nuclear-spin 
effects  contributing  to  the  hyperfine  structure  of  the 
6s’[  1/  2] ^  -6p’[3/2]2  xenon  ion  transition  produce  a 
total  of  18  isotopic  and  spin  split  components.  The 
hyperfine  splitting  constants  which  characterize  the 
variations  in  state  energies  are  only  known  for  a  limited 
number  of  energy  levels. 

It  is  not  presently  known  whether  hyperfine 
and  isotopic  constants  have  been  measured  for  the 
6^’[  1/  2]^  -6p’[3/2]2  transition.  Therefore,  it  is  not 
presently  possible  to  model  the  lineshape  to  determine 
temperature  and/or  pressure  broadening  parameters.  For 
velocity  measurements,  it  is  convenient  to  probe  opti¬ 
cally  accessible  transitions.  The  incomplete  knowledge 
of  the  isotopic  and  nuclear  spin  splitting  constants  does 
not  affect  the  determination  of  velocity.  A  convenient 
feature  of  this  line  is  its  proximity  to  the 
5J[3]7/2- 6p[2]5/2  ionic  transition  (~17  GHz  separa¬ 
tion).  This  fact  allows  for  the  same  optics  and  laser  to  be 
used  to  probe  either  transition  with  only  small  proce¬ 
dural  adjustments. 

Another  fortunate  feature  of  the 
6s’[  1/  2]j  -  6 p\  3/  2]2  transition  is  the  presence  of  a  rel¬ 
atively  strong  line  originating  from  the  same  upper  state 


(6i[3/2]j  -  6p'[  3/  2]2  transition  at  473.42  nm  [2,3]) 
which  allows  for  nonresonant  fluorescence  collection.  A 
nonresonant  fluorescence  scheme  is  preferred  where 
there  is  the  possibility  of  laser  scattering  from  surfaces. 

The  velocity  of  the  sample  volume  is  deter¬ 
mined  by  measurement  of  the  Doppler  shift  of  the 
absorbing  atoms.  If  an  absorber  has  a  velocity  compo¬ 
nent  u  along  the  axis  of  the  laser  beam,  it  will  absorb 
the  light  at  a  frequency  shifted  from  that  of  stationary 
absorbers.  The  magnitude  of  this  frequency  shift  5v12 
depends  on  the  speed  u  along  the  laser  beam  axis 


where  c  is  the  speed  of  light.  The  Doppler  shift  of  a  spe¬ 
cies  fluorescence  profile  away  from  the  line  center  v12 
of  stationary  absorbers  is  in  proportion  to  u  [4]. 

Apparatus 

The  measurements  in  this  effort  were  per¬ 
formed  in  Chamber  6  of  the  Air  Force  Research  Labora¬ 
tory  (AFRL)  Electric  Propulsion  Laboratory  at  Edwards 
AFB,  CA.  Chamber  6  is  a  nonmagnetic  stainless  steel 
chamber  which  has  a  1.8  m  diameter  and  3  m  length.  It 
has  a  measured  pumping  speed  of  32,000  1/s  on  xenon. 
Pumping  is  provided  by  four  single  stage  cryo-panels 
(single  stage  cold  heads  at  ~25  K)  and  one  50  cm  two 
stage  cryo-pump  (<12  K).  Chamber  pressure  during 
thruster  operation  is  approximately  5xl0"6  Torr,  cor¬ 
rected  for  xenon. 

Table  1 :  Nominal  Thruster  Operating  Conditions 


Anode  flow 

840  pg/s  (Xe) 

Cathode  flow 

98  pg/s  (Xe) 

Anode  potential 

250  V 

Anode  current 

0.820  A 

Keeper  current 

0.5  A 

Magnet  current 

0.75  A 

Heater  current 

3.0  A 

The  Hall  thruster  is  the  Busek  Company  BHT- 
200-X3  200  W  Hall  thruster  which  has  been  described 
elsewhere  [5].  Table  1  shows  the  nominal  operating  con¬ 
ditions  for  the  BHT-200-X3  thruster  for  this  test.  During 
thruster  operation,  the  parameters  shown  in  Table  1  are 
monitored  and  recorded  at  a  1  Hz  data  rate.  It  should  be 
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Thruster 


Fig.  1.  Side  view  diagram  of  thruster  within 
AFRL  chamber  6.  Also  shown  are  the  transla¬ 
tion  stages  and  the  laser  probe  beams.  Note  that 
the  fluorescence  collection  and  external  optics 
are  not  shown. 

noted  that  compared  to  previous  ion  velocity  measure¬ 
ments  [6],  the  magnetic  field  strength  was  reduced  due 
to  plasma  oscillations  which  affected  stable  thruster 
operation.  These  appear  to  be  a  function  of  thruster  age¬ 
ing  and  are  eliminated  by  a  slight  reduction  of  the  mag¬ 
netic  field. 

Figure  1  shows  a  side-view  diagram  of  the  Hall 
thruster  mounted  within  the  vacuum  chamber.  The 
thruster  is  mounted  on  a  three  axes  orthogonal  computer 
controlled  translation  system.  Figure  1  also  shows  the 
two  orthogonal  LIF  probe  beams  and  windows  through 
which  the  beams  enter  the  chamber.  Figure  2  shows  a 
top  view  of  the  laser  optical  train,  collection  optics,  and 
one  leg  of  the  external  probe  optics. 

The  laser  used  is  a  New  Focus  Vortex  tunable 
diode  laser.  It  is  capable  of  tuning  approximately  -50 
GHz  about  a  center  wavelength  of  834.7  nm.  The  10 
mW  beam  is  passed  through  a  Faraday  isolator  to  elimi¬ 
nate  feedback  to  the  laser.  The  laser  beam  then  passes 
through  several  beam  pick-offs  until  it  reaches  a  50-50 
beam  splitter  (BS)  where  it  is  split  into  two  beams  of 
equal  power.  The  first  beam,  the  axial  probe  beam 
shown  in  Figs.  1  and  2,  is  focused  by  a  lens  and  enters 
the  vacuum  chamber  through  a  window.  A  second  probe 
beam,  shown  in  Fig.  1  only,  is  directed  from  the  optical 
bench  via  a  periscope  apparatus  so  that  it  enters  the 
chamber  from  above  the  thruster  and  probes  the  velocity 


Fig.  2.  Top  view  diagram  of  the  laser  optical 
train  and  collection  optics.  Note  that  the  radial 
probe  beam  periscope  and  focusing  optics  are 
not  shown. 

perpendicular  to  the  first  probe  beam.  Each  probe  beam 
is  chopped  at  a  unique  frequency  by  choppers  Ch2  (2.8 
kHz)  and  Ch3  (2.0  kHz)  for  simultaneous  phase  sensi¬ 
tive  detection  of  the  two  fluorescence  signals. 

The  two  wedge  beam  pick-offs  (/BS)  shown  in 
Fig.  2  provide  portions  of  the  beam  for  diagnostic  pur¬ 
poses.  The  first  beam  pick-off  directs  a  beam  to  a  photo¬ 
diode  detector  (Dl)  used  to  provide  constant  power 
feedback  to  the  laser.  The  second  beam  is  divided  into 
two  equal  components  by  a  50-50  cube  beam  splitter. 
The  first  component  is  directed  to  an  Burleigh  WA-1500 
wavemeter  used  to  monitor  absolute  wavelength.  The 
second  component  is  sent  through  a  chopper  Chi  (1.3 
kHz)  and  through  a  low  pressure  xenon  hollow  cathode 
discharge  lamp.  The  lamp  provides  a  stationary  absorp¬ 
tion  reference  for  the  direct  determination  of  the  Dop¬ 
pler  shift  5v12.  The  second  pick-off  sends  a  beam  to  a 
300  MHz  free  spectral  range  Fabry-Perot  etalon  (F-P). 
This  instrument  provides  high  resolution  frequency 
monitoring  of  the  wavelength  interval  swept  during  a 
laser  scan. 
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The  fluorescence  collection  optics  are  also 
shown  in  Fig.  2.  The  fluorescence  is  collected  by  a  75 
mm  diameter,  300  mm  focal  length  lens  within  the 
chamber.  The  collimated  fluorescence  signal  is  directed 
through  a  window  in  the  chamber  side  wall  to  a  similar 
lens  that  focuses  the  collected  fluorescence  onto  the 
entrance  slit  of  125  mm  focal  length  monochrometer 
with  a  Hamamatsu  R928  photomultiplier  tube  (PMT) 
detector.  Due  to  the  1:1  magnification  of  the  collection 
optics,  the  spatial  resolution  of  the  measurements  is 
determined  by  the  geometry  of  the  entrance  slit  1  mm 
width  and  1.7  mm  height  as  well  as  the  sub-mm  diame¬ 
ter  of  the  probe  beams. 

The  laser  is  controlled  by  an  analog  ramp  sig¬ 
nal  generated  by  a  National  Instruments  E- series  data 
acquisition  board.  During  each  laser  scan,  the  data 
acquisition  card  records  the  absorption  and  two  fluores¬ 
cence  signals  using  three  lockin  amplifiers.  The  signal 
from  the  Fabry  Perot  etalon  photodiode  detector  (D3) 
signal  is  amplified  and  filtered  using  a  current  pre¬ 
amplifier.  The  output  of  which  is  also  recorded.  Typi¬ 
cally,  the  scans  span  6  GHz.  Each  scan  yields  four  traces 
of  several  thousand  points.  The  traces  are  then  stored  for 
analysis. 


The  data  analysis  procedure  consists  of  first 
normalizing  the  wavelength  interval  using  the  Fabry- 
Perot  etalon  trace.  Subsequently,  the  reference  absorp¬ 
tion  trace  is  rectified  by  subtraction  of  the  absorption 
baseline.  The  Doppler  shifts  of  the  LIF  signal  relative  to 
the  absorption  signal  are  then  calculated  by  two  meth¬ 
ods.  The  first  fits  an  estimated  lineshape  to  each  curve 
and  determines  the  Doppler  shift  from  the  calculated 
line  centers.  The  second  method  is  a  numerical  cross¬ 
correlation  which  compares  the  absorption  and  LIF  line 
shapes  to  determine  the  Doppler  shift.  It  is  much  less 
sensitive  to  signal  noise  and  is  more  mathematically  rig¬ 
orous  than  the  first  method.  However,  both  methods  typ¬ 
ically  produced  very  similar  results.  Measurement 
uncertainties  are  determined  from  the  differences 
between  the  two  methods  as  well  as  from  statistics  of 
repeated  measurements.  Typically  uncertainties  are 
greatest  near  the  main  discharge  due  to  low  signal 
strength  and  least  near  the  cathode  where  the  signal  is 
much  higher.  As  a  result,  the  LIF  scans  nearest  the 
anode  typically  integrated  for  20-25  min.,  while  the 
scans  taken  nearest  the  cathode  integrated  for  less  than  5 
minutes.  In  addition,  radial  measurements  near  the 
anode  discharge  exit  sometimes  exhibit  a  particularly 
wide  LIF  line  profile  which  cause  these  measurements 
to  have  the  greatest  variability.  For  the  measurements 
presented  here,  the  uncertainty  is  estimated  to  be 


Fig.  3.  Near  field  dimensions  of  the  BHT-200 
with  origin  of  the  coordinate  system  and  posi¬ 
tions  of  important  locations  noted. 

approximately  -50  m/s,  but  is  probably  lower  for  the 
axial  measurements. 

Figure  3  shows  the  near  field  geometry  of  the 
Busek  BHT-200  Hall  thruster.  The  locations  of  the  pro¬ 
truding  central  magnetic  pole  (nose  cone)  and  edges  of 
the  acceleration  channel  are  indicated  as  is  the  position 
of  the  cathode  exit.  The  cartesian  coordinate  system  and 
origin  used  in  these  measurements  are  also  shown  in 
Fig.  3.  The  coordinate  system  orientation  is  also  refer¬ 
enced  in  Figs.  1  and  2.  The  origin  is  at  the  tip  of  the  nose 
cone  due  to  the  ease  and  repeatability  with  which  this 
position  is  located.  All  measurements  are  identified 
using  these  coordinates.  Measurements  presented  in  this 
work  will  be  limited  to  the  Y-Z  plane.  Therefore  in  the 
anode  discharge  exit  region,  the  Y  axis  velocity  mea¬ 
surements  correspond  to  the  radial  velocities,  and  the  Z 
axis  measurements  to  the  axial. 

Experimental  Results 

Figure  4  shows  sample  LIF  and  rectified 
absorption  traces  used  to  determine  the  velocity  compo¬ 
nents  of  accelerated  xenon  neutrals  taken  at  the  cathode 
exit  at  Y  =  -33  mm  and  Z  =  +13  mm.  The  right-most 
lineshape  is  the  rectified  absorption  trace  (0  m/s),  the 
center  lineshape  is  the  axial  signal  (135  m/s),  and  the 
left-most  lineshape  is  the  radial  signal  (340  m/s).  The 
signal  at  the  cathode  exit  is  highest  signal  strength 
encountered.  Similar  signals  in  the  vicinity  of  the  anode 
discharge  exit  plane  required  lOx  longer  time  constants 
to  achieve  comparable  signal  to  noise  ratios  (SNR)  and 
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Fig.  4.  Typical  LIF  and  rectified  absorption 
trace  used  to  determine  velocity.  This  particular 
data  trace  is  from  the  cathode  exit  (Y  =  -33  mm, 

340  m/s;  Z  =  13  mm,  135  m/s). 

several  cases  required  30x.  This  limited  data  from  the 
anode  plume. 

Anode  Exit  Plane  Region 

Figure  5  shows  the  axial  and  radial  velocities  at 
the  anode  discharge  exit  plane.  The  upper  toroidal  por¬ 
tion  (see  Fig.  3)  is  the  only  portion  of  the  anode  plume 
where  velocity  measurements  are  available  to  the  geom¬ 
etry  of  the  protruding  central  magnetic  pole.  The  Xe  I 
velocity  is  nearly  completely  axial  along  the  outer 
radius  of  the  discharge  channel  (Y  =  16  mm).  The  axial 
velocity  component  diminishes  toward  the  central  mag¬ 
netic  pole  and  the  radial  velocity  component  appears  to 
increase. 

This  trend  is  mirrored  by  the  Gaussian  fit  pro¬ 
file  width  of  the  LIF  profile  shown  in  Fig.  6.  The  profile 
width  of  both  velocity  components  increases  consis¬ 
tently  with  proximity  to  the  protruding  boron  nitride 
capped  central  magnetic  pole.  Despite  lack  of  knowl- 
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Fig.  5.  Exit  plane  velocity  components  at  the 
exit  plane.  Note  error  bar  denoting  calculated 
uncertainty  in  the  velocity  measurements. 


edge  of  the  lineshape  function,  the  square  of  the  line- 
width  is  related  to  the  temperature  (e.g.  random  thermal 
energy)  [4].  It  therefore  appears  that  the  temperature  of 
the  neutrals  is  significantly  higher;  perhaps  as  much  as 
5x  across  the  width  of  the  acceleration  channel. 

Simply  using  the  radial  LIF  profile  width  to 
calculate  a  maximum  temperature  (simplistically 
assuming  purely  Doppler  broadening,  neglecting  all 
other  broadening  mechanisms),  the  temperature  ranges 
from  7,500  K  at  the  central  pole  to  1,500  K  at  the  outer 
radius.  As  a  reference,  the  lamp  discharge  profile  pro¬ 
duces  an  apparent  temperature  of  1,300  K.  These  esti¬ 
mates  ignore  the  complex  spectra  of  xenon  as  well  as 
any  pressure  broadening  mechanisms  and  represent  an 
absolute  maximum  temperature. 

It  appears  likely  that  neutrals  are  being  created 
by  ion  collision/neutralization  on  the  boron  nitride  sur¬ 
face  of  the  protruding  central  magnetic  pole.  These  neu¬ 
trals  can  not  be  differentiated  from  the  never  ionized 
neutrals.  The  combination  of  these  two  populations 
could  produce  broadened  velocity  distributions  similar 
to  those  seen  here.  Alternatively,  plasma  densities  may 
be  higher  near  the  central  pole  and  could  be  responsible 
at  least  some  of  the  broadening. 

Figure  7  shows  Xe  I  velocities  at  the  center  of 
the  acceleration  channel  (Y  =  12  mm)  extending  from 
the  exit  plane  (Z  =  -7  mm)  to  5  mm  beyond  the  tip  of  the 
central  magnetic  pole  (Z  =  +5  mm).  Both  the  axial  and 
radial  velocity  components  are  invariant  to  the  precision 
of  the  measurement.  Figure  8  shows  the  LIF  profile 
widths  of  the  axial  and  radial  LIF  signals  with  the  pro¬ 
file  width  of  the  absorption  trace  as  a  reference.  The 
radial  profile  width  decreases  as  would  be  expected 
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Fig.  7.  Velocity  components  along  the  accelera¬ 
tion  channel  centerline  (Y  =  12  mm)  extending 
from  the  exit  plane  (Z  =  -7  mm)  to  Z  =  5  mm. 

from  an  expanding,  cooling  flow.  Initially,  the  axial  pro¬ 
file  also  narrows,  although  only  slightly.  Beyond  the  tip 
of  the  central  magnetic  pole,  the  profile  increases  sud¬ 
denly.  It  does  not  appear  to  be  an  artifact  of  the  data,  and 
SNR  of  the  signal  is  not  significantly  different  in  this 
region.  A  possible  explanation  is  interaction  from  the 
neutral  flow  from  beneath  the  central  pole,  either  from 
the  opposite  portion  of  the  toroidal  channel,  or  from  the 
cathode  plume.  Alternatively,  the  effect  may  be  the 
result  of  a  gas  dynamic  shock  produced  by  the  interac¬ 
tion  of  an  over-expanded  anode  plume  interacting  with 
the  background  neutrals. 

Figure  9  contains  axial  and  radial  velocities  at 
the  outer  insulator  radius  (Y  =  16  mm)  extending  from 
the  exit  plane  (Z  =  -7  mm)  to  3  mm  beyond  the  central 
magnetic  pole  tip  (Z  =  3  mm).  The  data  is  very  similar  to 
that  of  the  acceleration  channel  center  line  data  shown 
in  Fig.  7.  The  axial  velocity  is  nearly  constant  at  approx- 
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Fig.  8.  LIF  profile  widths  (GHz)  for  both  axial 
and  radial  traces  along  the  centerline  of  the 
acceleration  channel  (Y  =  12  mm).  As  a  refer¬ 
ence,  the  profile  width  of  the  absorption  trace 
is  also  shown. 


Fig.  9.  Velocity  components  along  the  acceler¬ 
ation  channel  outer  radius  (Y  =  16  mm) 
extending  from  the  exit  plane  (Z  =  -7  mm)  to  Z 
=  3  mm. 
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Fig.  10.  LIF  profile  widths  (GHz)  for  both 
axial  and  radial  traces  along  the  outer  radius  of 
the  acceleration  channel  (Y  =  16  mm).  As  a 
reference,  the  profile  width  of  the  absorption 
trace  is  also  shown. 


imately  300  m/s.  The  radial  velocity  exhibits  more  vari¬ 
ation  than  the  previous  case,  but  the  behavior  is  similar 
with  an  average  velocity  of  approximately  75  m/s.  The 
corresponding  LIF  profile  widths  are  shown  in  Fig.  10. 
Similar  to  the  acceleration  channel  center  velocities 
shown  in  Fig.  8,  the  LIF  profile  width  of  the  radial  sig¬ 
nal  drops  with  distance  from  the  exit  plane,  typical  of  an 
expanding,  cooling  flow.  In  congruence  with  the  previ¬ 
ous  data,  the  axial  profile  width  is  near  constant  until  Z 
=  3  mm,  beyond  which  it  rises  by  nearly  25%. 


Cathode  Region 

The  hollow  cathode  consumes  approximately 
10%  of  the  total  thruster  propellant  flow.  The  cathode 
flow  supports  the  cathode  discharge  and  cools  the  cath¬ 
ode  to  prevent  evaporation  of  the  low  work  function 
electron  emitting  materials.  Due  to  the  restrictive  geom¬ 
etry  and  low  ionization  fraction,  the  neutral  xenon  den- 
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sity  is  significantly  higher  at  the  cathode  than  at  the 
main  discharge  exit,  by  an  estimated  20x. 

Figure  1 1  is  a  vector  plot  of  the  velocity  flow 
field  in  the  cathode  as  well  as  the  near  exit  plane 
regions.  At  the  cathode  exit,  the  neutral  velocity  is  380 
m/s.  The  remainder  of  the  flow  appears  to  be  a  gas 
dynamic  expansion.  The  peak  velocity  measured  in  the 
upper  right  hand  portion  of  the  cathode  flow  region  is 
approximately  610  m/s.  Interestingly,  the  maximum  is 
not  centered  on  the  initial  jet.  This  is  more  apparent  in 
Fig.  12  where  the  cathode  region  velocity  field  is  shown 
in  greater  detail.  The  individual  points  where  velocity 
data  is  available  are  denoted  by  the  arrow  origins.  The 
contours  show  the  interpolated  velocity  magnitude. 
Interestingly,  the  velocity  field  does  not  develop  sym¬ 
metrically  about  the  cathode  exit.  The  velocity  field  is 
skewed  away  from  the  anode  discharge.  The  cathode 
plume  appears  to  be  interacting  with  the  anode  plume. 

Some  evidence  of  the  interaction  of  the  anode 
and  cathode  plumes  is  shown  in  Figs.  13  and  14.  These 
show  the  Z  and  Y  lineshape  profile  widths.  (Note  in  pre¬ 
vious  discussion,  these  were  denoted  as  axial  and  radial, 
respectively.  This  terminology  in  the  cathode  plume  is 
not  appropriate  due  to  the  geometry.)  Figure  13  contours 
show  the  Z  axis  LIF  profile  width  as  it  varies  with  the 
cathode  plume  expansion.  It  shows  an  expansion  where 
the  profile  width  is  lessening  as  the  flow  expands  and 
cools.  There  appears  to  be  a  sudden  rise  downstream  of 
the  cathode;  however,  this  feature  is  due  to  a  single  data 
point  and  may  be  questionable. 

Figure  14  shows  the  Y  axis  LIF  profile  width  in 
the  near  cathode  region.  From  the  interpolated  profile 
contours,  it  appears  that  the  cathode  plume  contains  sev¬ 
eral  features,  possibly  shocks,  which  turn  the  flow  to  the 
right.  Taken  together,  Figs.  13  and  14  imply  that  the 
cathode  plume  is  interacting  with  either  the  anode 
plume  or  the  background  neutrals  in  the  vacuum  cham¬ 
ber  (background  pressure  ~5xl0"6  Torr).  This  interac¬ 
tion  appears  as  a  select  broadening  of  the  LIF  profiles.  It 
may  indicate  that  there  is  mixing  of  varied  velocity  dis¬ 
tributions,  or  that  a  gas  dynamic  shock  is  present.  Due  to 
the  sparseness  of  the  available  data,  there  is  consider¬ 
able  uncertainty  as  to  what  is  producing  the  increased 
LIF  profile  width,  but  it  indicates  that  the  anode  and 
cathode  neutral  plumes  exhibit  unexpected  complexity. 

Conclusions 

This  effort  demonstrates  neutral  xenon  velocity 
measurements  in  the  exit  plane  and  cathode  regions  of  a 
200  W  low  power  Hall  thruster.  The 


Fig.  11.  Neutral  xenon  vector  flow  field  in  the 
thruster  exit  and  cathode  regions. 


Z  (mm) 


Fig.  12.  Neutral  xenon  flow  field  in  the  cathode 
region.  Contours  denote  velocity  magnitude. 

6s’[l/2]j  -6p’[3/2]2  xenon  transition  used  in  this 
study  is  optically  accessible  using  relatively  inexpensive 
diode  lasers  and  is  approximately  17  GHz  from  the 
5d[3]7/2-6p[2]5/2 ionic  transition  previously  used  to 
measure  xenon  ion  velocities  [1,6].  Although  not  all 
transition  spectral  constants  are  known  for  either  transi- 
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Fig.  13.  Cathode  plume  Z  axis  LIF  profile 
width  contours  with  velocity  vectors  over- 
layed.  Note  the  general  cooling  trend. 


Z  (mm) 

Fig.  14.  Cathode  region  Y  axis  LIF  profile 
width  contours  with  velocity  vectors  over¬ 
lay  ed.  Note  the  apparent  interaction  with  the 
anode  plume 

tion,  their  proximity  allows  for  the  same  experimental 
apparatus  to  be  used  to  acquire  both  neutral  and  ionic 
velocity  data  with  minimal  adjustments. 

At  the  main  discharge  exit  plane,  there  is  a  pro¬ 
nounced  neutral  velocity  away  from  the  protruding  cen¬ 
tral  magnetic  pole.  The  axial  velocity  is  nearly  invariant 
(-300  m/s)  downstream  until  just  beyond  the  central 
magnetic  pole  where  a  sudden  increase  in  the  LIF  pro¬ 
file  width  occurs.  It  is  not  precisely  understood  what 
occurs  at  this  location,  but  it  is  possible  that  it  is  caused 
by  mixing  with  the  cathode  plume,  or  with  the  back¬ 
ground  neutrals.  It  is  also  possible  that  this  feature  is  the 


Fig.  15.  Photograph  of  the  BHT-200  Hall 
thruster  plume  operating  at  nominal  conditions 
with  various  features  identified. 


result  of  a  shock;  however,  there  is  insufficient  data  to 
confirm  any  of  these  possibilities. 

Due  to  the  low  degree  of  ionization,  the  cath¬ 
ode  plume  has  the  highest  concentration  of  xenon  neu¬ 
trals.  The  cathode  plume  exhibits  considerable 
acceleration.  The  velocity  measured  at  the  cathode  exit 
is  approximately  380  m/s  and  the  peak  velocity  mea¬ 
sured  further  downstream  is  610  m/s.  The  cathode 
plume  appears  to  be  slightly  skewed  by  interaction  with 
the  anode  plume.  Using  the  LIF  profile  widths  as  a  met¬ 
ric,  there  may  be  several  oblique  shocks  where  the  cath¬ 
ode  and  anode  plumes  interact.  However,  the  sparse  data 
indicates  that  caution  should  be  taken  in  this  analysis. 

Figure  15  is  a  photograph  of  the  BHT-200-X3 
Hall  thruster  at  AFRL.  Previous  ion  velocity  surveys 
have  determined  that  the  visibly  luminous  portions  of 
the  plasma  correspond  to  locations  with  mixing  of  sepa¬ 
rate  ion  velocity  populations  [6].  The  luminous  areas 
nearest  the  main  discharge  and  visible  cathode  jet  quali¬ 
tatively  appear  to  correlate  with  the  LIF  profile  width 
features  which  are  postulated  to  be  gas  dynamic  shocks. 
It  is  apparent  that  more  measurements  and  analysis  are 
required  before  these  artifacts  are  fully  understood. 

Cappelli  et  al.  performed  a  series  of  ground 
state  neutral  density  measurements  using  self  absorption 
in  the  vacuum  ultra  violet  [7].  These  measurements  sup¬ 
port  the  supposition  that  there  is  considerable  interac¬ 
tion  between  the  anode  and  cathode  plumes.  The  density 
measurements  also  show  an  anomalous  rise  in  neutral 
plume  density  directly  downstream  of  the  anode  cen- 


tered  near  the  cathode  plane.  This  feature  is  identified 
by  the  density  measurements  as  a  possible  gas  dynamic 
shock.  This  hypothesis  appears  to  support  the  conclu¬ 
sions  drawn  from  LIF  profile  widths  in  this  work.  It 
appears  obvious  that  the  neutral  flow  in  the  near  field  of 
a  Hall  thruster  is  a  complex  interaction  of  a  diffuse 
plume  form  the  main  discharge  and  a  more  concentrated 
jet  from  the  cathode.  The  degree  to  which  this  system 
interacts  with  the  ion  flow  from  the  main  discharge 
remains  an  open  question. 

Future  efforts  will  concentrate  on  understand¬ 
ing  how  to  interpret  the  differences  in  LIF  profile  width 
in  the  plume.  It  is  not  obvious  how  these  are  to  be  inter- 
pretated.  It  would  also  be  interesting  to  compare  these 
results  to  numerical  model  results.  Such  an  integrated 
study  of  the  near  anode  and  near  cathode  regions  could 
provide  valuable  insights  into  the  physics  of  the  flow 
interactions.  It  may  also  be  necessary  to  repeat  these 
measurements  for  a  transition  where  all  the  spectral  con¬ 
stants  are  known,  including  the  hyperfine  structure. 
Cedolin  pioneered  the  use  of  the  6s[3/2]2  -6p[3/2]2 
transition  at  823.2  nm  for  this  use  and  developed  a 
sophisticated  line  shape  model  [8]. 
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